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Abstract 



We have studied the relaxation and dephasing processes in a supercon- 
ducting quantum point contact induced by the interaction with an electro- 
magnetic environment. Based on a density matrix approach we obtain the 
rates for the dissipative dynamics as function of the transmission, the phase 
difference on the contact and the external impedance. Our calculation al- 
lows to determine the appropriate range of parameters for the observation of 
coherent oscillations in the current through the contact. 
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In a superconducting quantum-point contact (SQPC) the zero-voltage transport proper- 
ties are determined by the so-called Andreev states within the superconducting energy gap 
A, which are given by 



i?±(</.) = ±AVl-rsin2(0/2), (1) 

where (f) is the phase difference between the superconducting electrodes and r is the normal 
transmission coefficient for each conduction channel. Atomic contacts, produced by scanning 
tunneling microscope or the mechanically controllable break-junction (MCBJ) technique, 
have provided an almost ideal realization of an SQPC with a few conduction channels whose 
transmissions can be determined experimentally 0. Different methods for Andreev level 
spectroscopy have been suggested 0]. 

The Andreev states in a single channel SQPC constitute an interesting realization of 
a two-level system with well characterized parameters, which could be proposed as a solid 
state realization of a quantum qubit [Q. The phase difference through the contact can 
be fixed by the magnetic flux in a superconducting ring geometry like the one depicted 
in Fig. la. In addition, by varying the external flux one could prepare the system in a 
given linear combination of the Andreev states, in close resemblance to the case of a spin 
1/2 in a magnetic field. The phase would be, however, affected by fluctuations originated 
in the coupling of the ring with its electromagnetic environment, characterized by a finite 
impedance Z{uj). These fluctuations provide a source of relaxation and dephasing for the 
dynamics of our two- level system [Q. 

The aim of this work is to investigate this dissipative dynamics in order to determine 
the conditions for the observation of quantum- coherence effects in this system. 

For a single mode SQPC in the low bias regime one can neglect the continuous part of 
the quasi-particle spectrum and restrict the description to the subgap states, introducing 
the following Hamiltonian ^ 

i^o = A I cos(|) a, + r sin(^) a, I , (2) 



where are the Pauli matrix and r = Vl — r is the normal reflection amplitude of the 
contact. The Hamiltonian adopts the form (0) in the basis of ballistic states and becomes 
diagonal Hq = in the basis of the Andreev states. The ballistic and the Andreev 

states are represented in Fig. lb. 

In the ring geometry, the phase on the contact will be given by = 0o ~ ^4', where 
00 = i^e/h)^ (being $ the external magnetic flux) and S(j) represents the phase fluctuations 
due to the electromagnetic environment. If the impedance viewed by the contact is small 
compared to the resistance quantum Rq = h/Ae^ one can expand Eq.(^ to the flrst order 
in (50 to obtain the coupling between the subgap levels and the environment. Then, the 
Hamiltonian can be written as Ho{(j)) = Ho{4o) + He, where the coupling He = is 
proportional to the current operator /. Performing a rotation to the basis of the Andreev 
states, the current operator is given by 

i = — - I Vl — r (j^ — r sin(— ) cos(— ) o", I . (3) 
hE, 6.) r ^2^^2^'| ^ ^ 
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On the other hand, a generic environment can be represented by a set of L-C circuits, 
with phases (pn, such that 50 = J2n4'n- The relevant information of the environment is 
contained in the phase-phase correlation function [0] 

= y < 0„(t)0^(O) > = ^ f°° duj^^-^lcoth{^^)cosujt-ismu}t\ (4) 

RqJo UJ [ 2 J 

where Zt{uj) is the effective impedance of the circuit as seen from the contact and f3 is the 
inverse temperature. 

Projection operator techniques allow to derive a time-convolutionless generalized master 
equation (GME) ^ for the reduced density matrix p of the Andreev levels Up to second 
order in the coupling Hamiltonian this equation is given by 

m + ^[Ho,m] = -1 J^dr{[iJ{-r)p{t)] C{r) - [/, /)(t)/(-r)] C*(r)} , (5) 

where /(r) denotes the current operator (^) in the interaction representation. 

One can extract the evolution equations for the matrix elements of the reduced density 
projecting the GME (j^) in the {|+), |— )} basis of the Andreev states. Taking into account 

that P4.+ + P = 1, and defining pR and pi as the real and imaginary part of the off-diagonal 

matrix element we finally get the rate equations 

p+ = -{W+ + W')p+ + Apn , (6) 

PR = Vtpi + rip+ , (7) 
Pi = -{W^ + W- + r]A)pi - {Q + 2f],)p^ + 77^,(2^+ - 1) + T , (8) 



where p+ = p++, rj = rsin(</)o/2) cos(0o/2)/-\/l ~ ^ and = — EJ)/h. 

The coefficients entering in Eqs. (P) to (§) are the upward and downward transition 
rates which determine the system relaxation. The dephasing, i.e. the decay of the off 
diagonal elements, is also controlled by the coefficient A, while ^2,. plays the role of a level 
shift as discussed below. 

The explicit form of all these coefficients depends on the actual impedance Zt{uj). Typi- 
cally, this impedance would be determined by the measuring circuit, which could consist of 
a SQUID inductively coupled to the ring with the SQPC. We shall denote by M and L the 
mutual inductance between the ring and the SQUID and the self-inductance of the SQUID 
respectively. The SQUID is also characterized by a finite resistance R. The charging effects 
associated with the contact capacitance can be neglected in the low impedance regime that 
we are considering Under these conditions the real part of Zt will be given by 

If the environment has a large cuttoff frequency (7 = R/L ^ A//;,), we can make the 
Markov approximation in the evolution equations to (|^). In this case, the involved 
coefficients can be written as 
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T = «/„-!,-; (13) 



where 



_ 7rA SReZt(l]) 1-r 
° ~ 4 i?Q (1 -rsm2(0/2))3/2 ^ ^ 

is the downward transition rate at zero temperature, n{Q) is the Bose mean occupation 
number and Uk = 27r k/f3h are the Matsubara frequencies. 

The actual value of the coefficients is controlled by the environment parameters M, L 
and R. The experimental set up can be designed in order to fix these parameters in the 
appropriate range to observe quantum-coherence effects. Longer decoherence times would 
be obtained by reducing M as much as possible. However, SQUID parameters must satisfy 
certain constraints, imposed by the need to avoid hysteresis and thermal fiuctuations ||12|| . 
These considerations permit to estimate a minimal value of M of the order of 0.5nH, with 
L = O.lnH and R = 20fl Assuming that the superconducting material is aluminum 
{Aai ~ O.lSmel^), one obtains values of in the nsec~^ range. This rate fixes the order of 
magnitude of the coefficients ( p^Of ) to (p!3|). 

The system of differential equations is characterized by a set of eigenvalues which 
determine the typical times for the dissipative dynamics. In the limit of small impedance 
these eigenvalues are given by the following simple expressions 

Xi = -{w+ + W-), X2,3 = ~{w^ + w- + 2T]A)±i{n + nr). (i5) 

Notice that due to the presence of the environment, the Andreev levels are shifted ac- 
cording to E± ~ E± ±Qr/2. However, for the range of parameters we are considering this 
renormalization is small and can be neglected. 

Solving the equations (^-§1 at the lowest order in fi^^ (which is equivalent to implement 
the rotating wave approximation |T^) one realizes that the eigenvalues ( ]T5| ) are directly 



connected to the relaxation and dephasing rates by the relations Tj{ ~ — Ai and — 
-^e{\2,3}- 

The approximate evolution equations can then be written as 

P+ = ^(l-e-^-*)+P+(0)e-^^*, (16) 

i R 

PR = e"^°*(pi?(0) cos fit + pi{0) sin fit) , (17) 
Pj = e-^°\pi{0) cos fit - pr{0) sin fit) . (18) 

which clearly shows that the diagonal elements of the density matrix decay exponentially 
towards their equilibrium values on a time scale given by I/Tr while the non-diagonal 
elements perform oscillations with frequency fi which are damped on a typical time l/F^i. 
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This behavior is reflected in the evolution of the current mean value which is given by 



(/(t))=/o(l-2 




) + 2/oe-^«*( 




n) n 



e 



(pij(O) cos nt + p/(0) sin nt) , 



(19) 



where /q = e(A/^)^r sin is the equilibrium current at zero temperature. 

According to Eq. ([T9|) quantum coherence would manifest as damped Rabi oscillations 
in the mean current through the contact. These oscillations appear for an initial condition 
in which both the lower and the upper Andreev state are populated. An arbitrary initial 
condition can be in principle reached from the equilibrium situation by imposing a given 
time-dependent evolution to the phase. The amplitude of the oscillations are maximal for 
= TT where they become transmission independent. However, by varying the transmission, 
the frequency of the oscillations could be tuned in order to reach an experimentally accessible 
range. For instance, in the case of aluminum and transmission r = 0.99 the frequency fl 
would be of the order of 3 x 

The decay of these oscillations is controlled by the dephasing rate. Fig. 2 illustrates the 
dependence of the relaxation and dephasing rates on the phase 0o, the contact transmission r 
and the temperature. At zero temperature we have F^j = 2F£) = Wf^ , which means that both 
rates reach a maximum at 0o = ^ regardless of the contact transmission. This maximum 
decrease with transmission as a/I — r. At finite but low temperatures we have Tr = Wq{1 + 
and Yd = Wo{\ + 2r]'^kT/nQ). As expected both rates increase with temperature 
but the dephasing rate increases linearly instead of exponentially. The maximum of F;^ is 
located at (po = tt for all range of temperatures and behave as \/l — rcoth/5AVl — r. On 
the other hand. Yd develops a dip at 00 = tt and exhibits a double peaked structure located 



00 = 2arccos±Yl + (y (1 ~ ''")/2 — 1)/'*" for high temperatures. In this classical limit we 
get Yr = W,-2kT/nn and F^ = Yr^ + ^')- 



In summary, we have studied the dissipative dynamics of the Andreev states in a SQPC 
coupled to an electromagnetic environment characterized by a generic frequency-dependent 
impedance. Our results show that with the appropriate choice of the environment parameters 
one can reach values of the dephasing and relaxation rates within the nsec~^ range. A SQPC 
in a ring geometry thus provides an interesting realization of a solid-state two-level system 
in which quantum interence effects could be observed. 

We acknowledge useful discussions with C. Urbina, D. Esteve and A. Martm-Rodero. 
We are also indebted to F. Vazquez and A. Raguet for their stimulating comments. This 
work was supported by the Spanish CICyT under contract PB97-0044. 




5 



REFERENCES 



[1] A. Furusaki and M. Tsukada, Physica B 165&166, 967 (1990); C. W. J. Beenakker, 

Phys. Rev. Lett. 67, 3836, (1991). 
[2] E. Scheer et al, Phys. Rev. Lett. 78, 3535 (1997); E. Scheer et al, Nature 394, 154 

(1998). 

[3] L.Y. Gorelik et al. Phys. Rev. Lett. 81, 2538 (1998); A. Levy Yeyati, J.C. Cuevas and 
A. Martm-Rodero in "Photons and Local Probes", edited by O. Marti and R. MoUer, 
NATO ASI, Ser. E, Vol. 300, Kluwer (1995). 

[4] So far the most advanced results have been obtained on qubits based on flux states in 
small superconducting loops and on charge states in small superconducting islands. See 
Y. Nakamura et al., Physica B 280, 405 (2000) and V. Bouchiat et al., J. of Supercon- 
ductivity 12, 789 (1999). 

[5] Other sources of dephasing for Andreev states like Flicker noise have been studied 
recently by N.I. Lundin and Y.M. Galperin, Phys. Rev. 63, 094505 (2001). 

[6] D. A. Ivanov and M. V. Feigel'man, Phys. Rev. B 59, 8444 (1999). 

[7] G. L. Ingold and Yu. V. Nazarov, in Single Charge Tunneling, edited by H. Grabert and 
M. H. Devoret, NATO ASI, Ser. B, Vol. 294 (Plenum, New York, 1991). 

[8] T.M.Chang and J.L.Skinner, Physica A 193 (1993), 483. 

[9] Notice that our analysis is restricted to the subgap states neglecting coupling with the 
continuous part of the spectrum. Other interactions, such a electron-phonon coupling. 



would introduced a broadening of the subgap states |10]. For aluminum atomic contacts 



this broadening is, however, small [O] and would not introduce important corrections 

to the relaxation and dephasing rates. 
[10] A. Martin-Rodero, A. Levy Yeyati and F. J. Garcia- Vidal, Phys. Rev. B 53, 8891 (1996); 

D. Averin and H. T. Imam, Phys. Rev. Lett. 76, 3814 (1996). 
[11] M. F. Goffman et al, Phys. Rev. Lett. 85, 170 (2000). 
[12] M. Tinkam "Introduction to Superconductivity", McGraw-Hill (1996). 
[13] C. Urbina, private communication. 

[14] C. W. Gardiner, Quantum Noise (Springer- Verlag, Heildelberg, 1991). 



6 



FIGURES 




FIG. 1. (a) Schematic circuit: a superconducting quantum point-contact (SQPC) inserted in 
a superconducting ring threaded by a magnetic flux Z{io) denotes the effective impedance seen 
by the contact, (b) Andreev levels with phase dependent energies E± (full lines). The dotted and 
dashed lines indicate the ballistic states. 
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FIG. 2. Relaxation and dephasing rates Tji and F/j as function of the phase </> for transmission 
coefficients r= 0.9, 0.95 and 0.99 (from top to bottom) and two different values of the parameter 
/?A. The downward transition rate at zero temperature Wq coincides in the graphic scale with Tr 
for (3A = 100. 
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